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ABSTRACT 

A photometric redshift analysis of galaxies in the field of the wide-separation gravi- 
tational lens MG2016+112 reveals a population of 69 galaxies with photometric red- 
shifts consistent with being in a cluster at the redshift of the giant elliptical lensing 
galaxy z — 1.00. The A's-band luminosity function of the cluster galaxies is well rep- 
resented by the Schechter function with a characteristic magnitude K* — 18.90+ ; 57 
and faint-end slope a = —0.60^0 331 consistent with what is expected for a passively 
evolving population of galaxies formed at high redshift, Zf > 2. From the total Ks- 
band flux of the cluster galaxies and a dynamical estimate of the total mass of the 
cluster, the restframe ifs-band mass-to-light ratio of the cluster is derived to be 
M/Lks — 27^7 (M/Lk s )q, in agreement with the upper limit derived earlier 
from Chandra X-ray observations and the value derived locally in the Coma cluster. 
The cluster galaxies span a red sequence with a considerable scatter in the colour- 
magnitude diagrams, suggesting that they contain young stellar populations in ad- 
dition to the old populations of main-sequence stars that dominate the if s-band lu- 
minosity function. This is in agreement with spectroscopic observations which show 
that 5 out of the 6 galaxies in the field confirmed to be at the redshift of the lensing 
galaxy have emission lines. The projected spatial distribution of the cluster galaxies 
is filamentary-like rather than centrally concentrated around the lensing galaxy, and 
show no apparent luminosity segregation. A handful of the cluster galaxies show evi- 
dence of merging/interaction. The results presented in this paper suggest that a young 
cluster of galaxies is assembling around MG2016+112. 

Key words: galaxies: clusters: individual: MG2016+112 - galaxies: elliptical and 
lenticular, CD - galaxies: evolution - galaxies: formation -galaxies: luminosity function, 
mass function - cosmology: observations 



1 INTRODUCTION 

1.1 Galaxy evolution in clusters 

The predominance of early-type galaxies in clusters com- 
pared with the field is a clear sign that the evolution of 
galaxies depends on environment, but the physical mecha- 
nism at work, in particular the question of when and how 
the most massive galaxies, the giant ellipticals, formed is a 
topic of ongoing theoretical and observational research. 

The traditional "monolithic" e lliptical galaxy formation 
scenario proposed bv lEggen et alJ lll962tl postulates a single 
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burst of star formation at high redshift, followed by passive 
stellar evolution. 

Early type galaxies in nearby and intermediate red- 
shift clus ters are observed to form a tight colour-magnitude 
relation feower et a,ljll992j: |Aragon-Salamanc a et aljll993l : 
iGladders et alJ Il99ct [Stanford et al.1 ll998Tl . or "red se- 
quence" which when combined with the close correlation 
between galaxy mass and m etallicity implied by the Mg2-er 
relation feender et aljll993) is explained naturally by a sin- 
gle, early episode of star formation [zf > 2). If more massive 
galaxies are more efficient at retaining supernova ejecta, they 
will have higher metallicities, and therefore redder colours 
jArimoto fc Yoshiilll987l : iKodama fc Arimotd|l997t) . 

There is however many pieces of evidence suggest- 
ing that the evolution of the entire cluster galaxy popu- 
lation is not as simple as expected from the monolithic 
galaxy formation scenario. The fraction of blue galaxies in 
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clust e rs inc r ease rapidly with redshift feutcher fc Oemlerl 
ll97St 11984 iRakos fc Schomberd Il995l) . High resolution 
imaging and spectroscopic studies show that a signifi- 
cant fraction of intermediate redshift blue galaxies are 
late-type s pirals with emission line s and/or post-starburs t 
signatures dPressler fc Gunnl l|l992h : lDressier et alJ (1999); 
iPoggianti et alJ dl999f) and that the cluster SO popula- 
tion is far less abundant at z ~ 0.5 than today. It has 
been proposed that these blue star-f orming galaxies are 
the progenitors of cluster SO g alaxies iDressler et alj|l997t 
Ivan Dokkum et alJll99gt l2000l) . In agreement with this ex- 
planation of the Butcher-Oemler effect, some SOs in the 
present epoch appear to con tain younger stellar populations 
ilKuntschner fc Davieslll99Sl) . 

Hierarchical structure formation models present a very 
different view of galaxy evolution in which galaxies assem- 
ble by merging of smaller stellar systems over a wide redshift 
range. In this scenario, the evolution of elliptical galaxies are 
governed by the time at which the bulk of the stellar popula- 
tions are formed, the era when the majority of mergers took 
place, and the amount of new star formation induced during 
each merger event. If the formation of the stars in an ellip- 
tical galaxy takes place over a broad cosmic time interval, 
then ellipticals at any redshift should exhibit a wide range of 
mass-weighted stellar ages. If however most stars in present 
day cluster ellipticals were formed in smaller disks at high 
redshift z > 1 and if little additional star formation took 
place in subsequent mergers, then the end product ellipti- 
cals could appear to be old and approximately coeval even if 
the bulk of m erging took place relatively late . The hierarchi- 
cal models of iKauffmann fc Charlotl <U99 8a ,v ) reproduce the 
colour-magnitude relation of ellipticals because more mas- 
sive galaxies are formed from merging of systematically more 
massive progenitors, which retain more metals when forming 
stars. 

Discriminating between the two formation scenarios de- 
scribed above is difficult, even when studying galaxies in 
cluste rs at high redshift. An increasing number of observa- 
tions {Stanford et al . 1997; lEUis et alJ 119971: IStanford et alJ 
Il998l: iRosati et al.lll999l : iNakata et al]l200ll) show that the 
evolution of the "red sequence" to redshifts z ~ 1 is consis- 
tent with a high formation redshift and a subsequent passive 
evolution favored by monolithic collapse, but none of them 
can exclude the possibility that cluster ellipticals formed 
from mergers of smaller galaxies as long as the bulk of star 
formation took place at much larger redshifts and there was 
little star formation in the subsequent merging process. 

One way to distinguish between the two scenarios is by 
studying the evolution of the mass function of the cluster 
galaxies. Hierarchical galaxy formation models predicts the 
mass assembly of galaxies to take place over a long time 
scale (z <C 2) and to be quite dec oupled from star forma- 
tion ijKauffmann fc Charlodll998a) . If this picture is true, 
a strong evolution of the stellar mass function with redshift 
is expected. This evolution is best studied in the .Ks-band 
because it is relatively unaffected by ongoing/recent star- 
forma tion and hence directly measu res the growth of galaxy 
mass (IKauffmann fc Charlotlll998bll . 

Locally, the cluster galaxy luminosity function have 
been studied in great de tail both at optic al and Near In- 
frared (NIR) wavelengths. lGoto et al.l <|2002T) derive the com- 
posite luminosity function of 204 "Sloan Digital Sky Survey 



Cut & Enhance Galaxy Cluster Catalog" in the redshift 
range range z = 0.02 to z = 0.25 in the five SDSS bands 
U* , g* , r*, i* and z* , and find that the faint-end slope of 
the luminosity function becomes flatter toward the redder 
wavebands, consistent with the hypothesis that the clus- 
ter luminosity function has two distinct underlying popula- 
tions; a population of bright ellipticals with a gaussian-like 
luminosity distribution that dominate the bright end, and a 
population of faint blue star-forming galaxies with a steep 
power-la w like luminosity distrib ution, that dominate the 
faint-end. Ide Propris et al] Jl99St) derive the NIR (H-band) 
luminosity function of the Coma cluster galaxies, and find 
results consistent with the above picture of a population of 
bright red galaxies and large population of faint blue dwarf 
galaxies. 

HST imaging of high redshift cluster galaxies have re- 
vealed a pop ulation of luminous "red" merg ers in high red- 
shift clusters llvan Dokkum et all200oLl200 ll) . The existence 
of these galaxies which follow the same colour-magnitude 
relation as the cluster ellipticals (but with a slightly larger 
scatter) supports the hierarchical galaxy formation scenario. 

The bright (massive) end of the luminosity function is 
expected to steepen and shift to fainter magnitudes at high 
redshift, as the old bright cluster galaxies begin to break up 
into building blocks. 

In contrast to the bright cluster galaxies, the faint clus- 
ter galaxies seem t o have a much greater d iversity of star 
formation histories dKodama fc Boweil feoOl). The Butcher- 
Oemler effect in intermediate redshift is an example of this. 
If this progression of activity to lower-mass galaxies as the 
universe ages is a consequence of the faint galaxies forming 
at later cosmic times than their massive counterparts, then 
distant clusters should have a luminosity function with a de- 
clining faint-end slope in contrast to r ising faint-end slope 
of local clusters llde Propris et al1 IT998l, 

The evolution of the if-band luminosity function in 
the redshift range z = 0.1-0.9 has been studied by 
Propris et al] a!999f) . The evolution of the characteris- 
tic magnitude K* of the galaxies, is found to be consistent 
with what is expected for a passively evolving population of 
galaxies formed at Zf > 2. The data is not deep enough to 
constrain the faint-end slope which is fixed at a = —0.9 (the 
value derived in the H-band for the Coma cluster). The au- 
thors find no significant dependence on the X-ray luminosity 
of the clusters. 

Naka ta~et al] feOOll) extend the study to z = 1.2 by de- 
riving the K band luminosity function for the cluster around 
the radio galaxy 3C324. These authors also choose to fix the 
faint-end slope at a = —0.9 and find consistency with what 
is expected for passively evolving stellar populations formed 
at Zf > 2. 

The suggested passive evolution of K* with redshift to 
z = 1.2 is a nice confirmation of the passively evolving stellar 
population found from other methods (such as the colour- 
magnitude relation) , but does not help discriminate between 
the monolithic and hierarchical formation scenarios, since no 
steepening of the bright end of the luminosity function (a 
deficit of the brightest galaxies) has been observed and no 
attempt has been made to measure the faint-end slope. 
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1.2 Previous observations of the MG2016+112 

system 

MG2016+112 is a wide separation gravitational lens first 
identified by lLawrence et alJ (Il984l) to be comprised of 
at least three lensed images of an active galactic nucleus 
(AGN) at redshift z — 3.269. The primary lensing object 
is a giant ell iptical galaxy a t z = 1.004 (hereafter "th e 
lens galaxy". ISchneider et alJ (Il986l) : ISoucail et al.1 fcOOlt) : 
iKoopmans fc Treul i2002Tl 'l. Wide separation gravitational 
lens systems have in many cases been found to contain 
galaxy clusters or groups that contribute significantly to the 
large separation angles. A nice property of clusters discov- 
ered in this way is that they are mass selected, and there- 
fore not biased toward the old relaxed systems that clusters 
selected in X-rays or in the optical/NIR are (since their se- 
lection depends on the presence of thermalized gas and/or 
central concentrations of luminous old red galaxies). 

lHattori et al. (1997) reported the detection of resolved 
X-ray emission at z ~ 1 centered on the lens galaxy. This re- 
sult was based on the detection of an X-ray line at 3.35 keV 
in the ASCA spectra of the system, which was interpreted as 
Fe K emission originating from a cluster at z ~ 1, and on the 
detection of "diffuse" X-ray emission near the giant elliptical 
in ROSAT HRI observations. Infrared and optical observa- 
tions of the field rev ealed no excess of galaxies around the 
giant elliptical lens Schneider et aljll985t lLangston et alJ 
1 19911: lLawrence et alill993l) .' This led to the suggestion of 
a "dark cluster" associated with the lens, in which a clus- 
ter sized dark matter and hot gas overdensity exists with 
few optically bright galaxies (M/L v > 2OOOM /L ). This 
conclusion was puzzling as the observed iron line in the X- 
ray spectra indicated a near solar metallicity for the cluster, 
which implies a long history of star-formation in the region. 
iBenitez et al.l (|l999j) re-analyzed the ROSAT HRI data and 
found a faint elliptical X-ray source which was coincident 
with a red galaxy overdensity. Deep V ,1 and K band Keck 
imaging revealed a red sequence of galaxies, with colours 
consistent with being at redshift z ~ 1. 

ISoucail et al.l ll200 if) performed a spectroscopic survey 
of 44 galaxies in the field and found an excess of 6 red galax- 
ies securely identified to be at z ~ 1, though not very cen- 
trally concentrated around the giant elliptical galaxy. For 
these six galaxies a mean redshift of z = 1.005, a velocity 
dispersion of a = 77lt 4 gokms _1 , a dynamical mass esti- 
mate of 2.8 x 10 14 /i5oMq and a V-band mass-to-light ratio 
of M/L v = 215t^ 8 /i5oM /L Q was derived. 

It is interesting to note that a high fraction of the z = 1 
galaxies (5 out of 6) have emission lines in their spectra, 
indicating the presence of young stellar populations. This is 
a much larger fraction than in low redshift clusters which are 
dominated by old quiescent galaxies with no emission lines, 
and also much higher than in other high redshift (z > 1) 
clusters, wher e the fraction of emission line galaxies is found 
to b e 10-15% llStanford et al.lll997l : Hlosati et al.lll999ft . 

IClowe et alJ (1200 ll) performed a weak lensing analysis of 
the field and detected a mass concentration 64" northwest 
of the lensing galaxy, and were able to rule out (at the 2a 
level) a mass pea k centered on t he len sing galaxy of the size 
suggested by the lHattori et al.l dl997h X-ray observations. 

IChartas et al.l (1200 it) recently acquired a 7.7 ksec ob- 
servation of the system with Chandra. These observations 



can account for all the X-ray emission as originating from 
the lensed images and additional X-ray point sources in the 
field, and place an 3a upper limit on the 2-10 keV flux and 
luminosity of diffuse cluster emission of 1.6 x 10~ 14 ergs 
s _1 cm -2 and 1.7 x 10 44 ergs s _1 , respectively, and an es- 
timated upper limit on the mass M 500 < 4.5 x W 14 h 50 M Q 
and mass-to- light ratio of M/Lv < 190Mq /Lq. 

Rather than being an X-ray luminous "dark-cluster", 
MG2016+112 thus turned out to be a spectroscopically con- 
firmed cluster of galaxies, detected neither in X-rays nor in 
weak lensing maps. 

In this paper we revisit the system through deep NIR 
and optical observations. We identify 69 galaxies in the 
field with photometric redshifts consistent with being cluster 
galaxies at redshift z — 1 and derive their colour-magnitude 
sequences and their luminosity function. 

The outline of the paper is as follows: In section [5] we 
describe the observations and data reduction. In section |3] 
we describe the construction of a multi colour photometric 
catalog of galaxies in the field (object detection, matching, 
photometry, star/galaxy separation). In section0]we discuss 
the photometric redshift analysis. The results are presented 
in section and summarized in Section |S| 

Throughout this paper we assume a flat = 0.3, Qa = 
0.7, ho = 0.70 cosmology unless otherwise noted. 



2 OBSERVATIONS AND DATA REDUCTION 

The observational basis of this work consists of very deep 
NIR J and Ks band observations carried out with the 
ISAAC instrument (proposal 63.O-0379) in service mode 
between May and July 1999, and moderately deep opti- 
cal B, V, R, and I band observations carried out with the 
FORS1 instrument (proposal 65.O-0666) in service mode in 
July 2000. Both FORS1 and ISAAC are mounted on the 
VLT/ANTU telescope. ISAAC employs a 1024 2 Hawaii ar- 
ray, with a pixel scale of 0'.'147 and a field of view (FOV) of 
2!5x2^5. FORS1 employs a 2048 2 CCD (TK2048EB4) which 
with the high resolution collimator used for the present ob- 
servations has a pixel scale of O'.'l and a FOV of 3.'4 x 3f 4. 

The data was originally obtained for other purposes 
than studying the cluster galaxy population. The NIR 
dataset was designed to perform a weak lensing analysis on 
the background galaxies. The optical dataset was designed 
to study the central bright gravitationally lensed QSO and 
is therefore not deep enough to detect the fainter cluster 
galaxies. The derived magnitudes (or upper limits) however 
serve as powerful input for the photometric redshift analysis 
which separates the cluster galaxies from the foreground and 
background field galaxy population. All observations were 
performed in excellent seeing conditions. The observing log 
is given in Tab.Q 

2.1 Reduction and calibration of NIR data 

The deep near-infrared imaging consists of 70 images of 2 
minutes integration time (2 hours 20 minutes) obtained in J 
in 3 independent sequences, and 275 images of 1 minute each 
(4 hours 35 minutes) obtained in Ks in 5 sequences. Some 
of them were taken in photometric conditions while others 
not, so subsequent re-scaling to the photometric ones was 
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Table 1. Log of the VLT/UT1 observations 



Instrument 


Filter 


NExp 


Exp 


Seeing 


Completeness 


ISAAC 


Ks 


275 


60 


0'.'4 


22™ 


ISAAC 


J 


70 


120 


0"4 


23 m 


FORS1 


I 


2 


200 


0'.'4 


24 m 


FORS1 


R 


2 


150 


C'4 


24™ 5 


FORS1 


V 


1 


180 


or.' 4 


26 m 


FORS1 


B 


5 


300 


or.' 5 


27^0 



NExp: number of exposures, Exp: exposure time in seconds 



applied in the final combination. Data red uction was per - 
formed with the "eclipse" software package (Dcvillard 1997) 
and IRAF routines. Flat-fields were produced from series of 
twilight sky fiats. To take into account the rapidly varying 
NIR background the sky value in each pixel was calculated 
as a running median of the pixel value in the exposures taken 
immediately before and after a given exposure. 

To compensate for the dithering pattern and optical 
distortion the individual frames were transformed to a ref- 
erence frame using a 3rd order geometrical transformation 
calculated from the position of 30-40 stars in each frame. 
Finally, the images were averaged with some clipping to 
remove residuals defects. Photometric calibration was ob- 
tained through the observations of standard calibration stars 
during the photometric nights of the programme. 



3 OBJECT DETECTION, PHOTOMETRY 
AND MATCHING 

Object detection was done in the Ks ban d, using the SEx- 
tractor software feertin fc Arnoutslll99r3) with default pa- 
rameter settings. We searched for the detected objects in 
the remaining frames using the SExtractor ASSOC option, 
which assumes that the frames being matched have been 
registered to the same pixel coordinate system. 

The transformations of the combined images to the Ks- 
band image were derived using the IRAF tasks geoform. 

Transforming the optical images to the Ks-band im- 
age includes a rebinning (due to the different pixel scales of 
the two instruments) in addition to the distortion transfor- 
mation. To avoid this image degradation, as well as other 
potential problems associated with non-linear image trans- 
formations, such as blurring of persistent bad pixels and cos- 
mic ray events, corruption of the background statistics etc. 
we transformed the Ks-hand object list xy-pixel coordinates 
to match the reference frames of the other images using the 
IRAF task geoxytran, rather than transforming the images 
themselves. In this way the photometry is performed on the 
original, untransformed images. 

For the total Ks band magnitudes of the objects we 
used SExtractors maq_ best which has bee n shown to be ac- 
curate to about 0^05 JSmail et alJl200lTl and colours were 
derived in fixed apertures of approximately 12 kpc at z — 1 
(10 pixels in the ISAAC images, 15 pixels in the FORS1 
images) . 



2.2 Reduction and calibration of optical data 

The individual optical data frames were bias subtracted and 
flat fielded using the standard ESO pipeline a nd cleaned for 
cosm ic rays using Laplacian edge detection (Ivan Dokkuml 
l200ll) . In some of the frames the background was found to 
have large scale gradients over the field. To properly ac- 
count for this we ran all the fram es through the SExtractor 
software feertin fc Arnoutsil99tJl with options set to save a 
full resolution interpolated background map, which was then 
subtracted from the science frames (large scale background 
variation are well modeled by SExtractor) . Cosmetic defects 
in areas with no objects were patched by interpolation of the 
background in surrounding annuli. This was necessary since 
some wavebands had too few exposures for deviant pixel 
rejection. Finally the individual frames were geometrically 
transformed to the same reference frame and combined using 
the IRAF tasks geoform, geotran and imcombine. 

The combined frames were calibrated using the atmo- 
spheric extinction and color corrected zero-points quoted 
on the ESO webpage, and corrected for Galactic extinction 
us ing the parameterization of the G alactic extinction law 
of ICardelli. Clayton fc Mathisl ( 1989) , assuming the mean 
Galactic value for the total-to-selective extinction in the 
V band R v = {Rv) MW = 3.1 and E(B- V)=0.22 (from 
interp olation of DIRBE/IRAS dust maps, ISchleeel et alJ 
1998). Finally, all the f rames were WCS cal ibrated using 
the USNO-A2.0 catalog llMonet fc et alJll998ft . 



3.1 Completeness 

To estimate the completeness magnitude of the observations 
we ran SExtractor object detection independently on all the 
frames. In Fig.^we plot the number of detected objects as a 
function of Ks and J band magnitude. From visual inspec- 
tion of the turnover magnitudes of the distributions in Fig.0 
we estimate that the Ks band observations are complete to 
about 22™0 and the J band observations are complete to 
about 23 m 0. The estimated completeness magnitudes of all 
the wavebands, including the optical are listed in Tab.0 



3.2 Star/galaxy separation 

SExtractor assigns to each object a parameter class_star 
which ranges from to 1. When the FWHM is well deter- 
mined stars obtain values close to 1, while extended objects 
obtain values close to 0. A main concern is contamination 
of the galaxy sample by stars, especially since some of them 
(e.g. M stars) would have the same colours as early-type 
galaxies at the cluster redshift. We therefore adopt a very 
conservative threshold in the classstar parameter of 0.05. 
Objects with class_star larger than this value in the Ks band 
or J band are labeled stars. Now we can construct a multi 
colour catalogue of galaxies consisting of objects which (i) 
have Ks < 22 (ii) are detected in at least two bands and 
(iii) are classified as a galaxy by means of the classstar pa- 
rameter. 
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Figure 1. Binned representation (binsize is 0™5) of the number 
of detected objects as a function of magnitude in the Ks (full 
curve) and J band (dashed curve). The data are estimated to be 
complete to Ks=22 ra and J=23™0 
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Figure 2. Photometric redshifts distribution of the galaxies in 
the field of MG2016+112 derived using hyperz. 



4 PHOTOMETRIC REDSHIFT ANALYSIS 

To separate the cluster galaxies from the foreground and 
background galaxies we estimated their redshift s using the 
public hyperz photometric redshift code of iBolzonella et all 
(1996). The code uses the GISSEL98 (Galaxy Isochrone Syn- 
t hesis Spectral Evolut i on Lib rary) spectral evolution library 
of lBruzual fc Charlou il993T) to build synthetic spectra with 
7 different star formation histories which match the observed 
properties of local galaxies from E to Irr type: a constant 
star-forming rate and six exponentially decaying star for- 
mation rates with time scales from 1-30 Gyr. The models 
assume solar metallicity and a Millar- Scalo IMF, and inter - 
nal reddening is considered using the lCalzetti et ail j2000f) 
model with Av varying between 0™0 and 1™2. Objects not 
detected in a waveband are assigned zero flux (with an er- 
ror corresponding to the limiting magnitude) rather than 
just ignored in the fit. This "upper limit" is important for 
the photometric redshift determination and is the reason 
for including the rather "shallow" optical data in the anal- 
ysis. Fig. |21 shows a histogram of the derived photometric 
redshifts of all the galaxies in the field. There is a distinct 
peak at redshift Zphat ~ 1, and quite a few galaxies with 
Zphot < 0.2. 

In Fig. |3] we plot the derived photometric redshifts of 
all the galaxies in the field versus their Ks-h&nd magnitude. 
The open squares mark the galaxies which have been spec- 
troscopically confirmed to be at redshift z = 1 (hereafter 
spectroscopic members. ISoucail et al.l {2001)). The full line 
marks the redshift of the lensing galaxy at z = 1 while the 
dashed lines are the limits z p hot = 0.7 and z v y tot = 1.2 that 
encompasses all the galaxies in the peak around z v hot = 1> 
including the 5 spectroscopically confirmed ones. Fig. |S| re- 
veals that most of the z p hot J; 0.2 are faint Ks > 20 galax- 
ies which may have been misclassified as dusty low redshift 
galaxies due to poor photometry. 




Figure 3. Photometric redshifts of all galaxies in the field of 
MG2016+122 versus their Ks band magnitudes. The error bars 
are la from hyper-z. The squares mark galaxies wh ich have been 
spect roscopically confirmed to be at redshift 2 = 1 ISoucail et alj 



We classify the 69 galaxies within the limits 0.7 < 
Zphot < 1.2 as probable cluster galaxies (hereafter photo- 
metric members). The choice of limits is a trade off between 
maximum cluster galaxy completeness and minimum con- 
tamination from field galaxies. We choose the lower limit to 
be Zphot > 0.7 rather than e.g. z p hot > 0.8 to avoid excluding 
spectroscopic members, and other cluster galaxies with simi- 
lar colours. Young stellar populations in the cluster galaxies 
can influence their broadband colours and derived photo- 
metric redshifts due to the discrete nature of the spectral 
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Figure 4. Composite Ks, Js, /-band image of the field of MG2016+112. The field is 2.5'x2.5'. The pixel scale is 0"l47pix~ 1 . The 
circles mark the 69 "photometric members" of the cluster at z = 1. The squares mark the "spectroscopic members" of the cluster. North 
is up, East is to the left. 



library used by hyperz. Since a large fraction of the spectro- 
scopically confirmed cluster galaxies show evidence of recent 
star formation it is likely that this is also the case for other 
cluster galaxies. Some galaxies in the 0.7 < z p hot < 0.8 
interval may be truly foreground, but many could also be 
cluster galaxies at z = 1 with recent star formation. To 
test the stability of the photometric redshifts derived with 
hyperz versus the photometric errors we did the following: 
For each Am=0.5 mag bins in the observed range of mag- 
nitudes Ks = 16 — 22, we generated a catalogue of 1000 



galaxies at z — 1 for the same bandpasses and limiting mag- 
nitudes as for the data of MG2016+112. The catalogue was 
generated using the make_catalog code in the hyperz pack- 
age. The galaxies are randomly drawn from the 7 template 
spectral types (E to Im) used for estimating the photometric 
redshifts. We then apply hyperz to the catalog to estimate 
the photometric redshifts. In Fig. we plot the fraction of 
the input [z = 1) galaxies that have derived photometric 
redshifts in the range 0.7 < z v hot < 1-2 as a function of 
ifs-band magnitude. From this figure it can be seen that 
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Figure 5. Simulated fraction of redshift 2 = 1 galaxies with 
derived photometric redshifts in the range 0.7 < z p i lot < 1.2 as a 
function of Ks band magnitude 

hyperz recovers 100% of the galaxies down to magnitudes 
of Ks = 18.5. At Ks — 20 it recovers about 90% and at 
Ks = 22 the recovery rate is down to 40%. 

The redshift evolution of the field galaxy luminosity 
function, and its cosmic variance is not known with suffi- 
cient accuracy to correctly incorporate the effects of pollu- 
tion from field galaxies in the completeness function analy- 
sis. Instead we will return to the issue of field contamination 
in Sec. l5.1l where we discuss the derived projected spatial dis- 
tribution of the cluster galaxies, and in Sec. 15.41 where we 
discuss their luminosity function. 



5 RESULTS 

5.1 Distribution of the cluster galaxies 

The distribution of the cluster galaxies (photometric mem- 
bers) on the sky is not very centrally concentrated. This is 
illustrated in Fig.0]where the overlay circles mark the photo- 
metric cluster members, and in Fig.^Jwhere the overlay rep- 
resents contours of the smoothed density distribution of the 
photometric member galaxies. The contours are 3-30u above 
the background density which is taken to be the density of 
galaxies in the HDFN with K 3 < 22 and 0.7 < z phot < 1.2 
l)Fernandez-Soto et alJlT999l ). The apparent morphology of 
the cluster is an elongated filamentary-like distribution sim- 
ilar to what i s observed in other z > 1 clusters (e.g . the 
Lynx clusters ijStanford et al]fl99 j]Jlj^ati_^^ and 
the RDCS J0910+5422 cluster l|Stanford et al.ll2002l n and 
consistent with what is expected for a "dynamically young" 
cluster in the process of assembly. The high significance of 
the overdensity structure, compared to the HDFN indicates 
that most of the photometric members are z = 1 cluster 
galaxies rather than random field galaxies that would be 
distributed more homogeneously across the field. Fig. il- 
lustrates that there is no apparent luminosity segregation 




Figure 6. Contours of the smoothed density distribution of 
photometric member galaxies overlayed the _R"s-band image of 
MG2016+112. The contours are 3-30<r above the background den- 
sity which is taken to be the de nsity of galaxies in the HDFN with 
K 3 < 22 and 0.7 < z phot < 1.2 <Fernandez-Soto et al.ll99flt) . The 
largest concentration of photometric member galaxies is centered 
on the central lensing galaxy, and is part of a "loose" elongated 
filamentary like structure. The pixel scale is 0Vl47pix~ 1 . North 
is up, East is to the left. 

of the cluster galaxies (i.e the brighter cluster galaxies are 
not more centrally clustered than the fainter galaxies) . This 
is consistent with the picture of a young cluster where the 
smaller galaxies have not yet had time to merge and trans- 
form into the population of massive galaxies found in the 
center of relaxed clusters. 

5.2 Properties of the cluster galaxies 

A full morphological analysis of the cluster galaxies requires 
higher resolution imaging data, and is beyond the scope of 
this paper, but qualitative inspection of Fig.0]sliows a wide 
variety of galaxies types: elliptical, spirals and even three 
or four mergers. Merging cluster galaxies are very rare in 
the cores of local relaxed clusters. The significant fraction 
of merging cluste r galaxies found here, and in o ther high 
redshift clusters (Ivan Dokkum et al J 12000. 2001) is direct 
evidence of the hierachical formation scenario. 

5.3 Colour-magnitude relation 

In Fig.|H|and Fig.|5|we plot the J — Ks versus Ks and R—Ks 
vs Ks colour-magnitude relations of galaxies in the field of 
MG2016+112. Filled symbols mark galaxies within 0.5 Mpc 
of the lensing galaxy. We label this region the "inner region" . 
Open symbols mark galaxies farther away than 0.5 Mpc from 
the lensing galaxy. We label this region the "outer region". 
The area of the outer region is approximately the same as the 
inner region. Also indicated are predicted colour-magnitude 
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Figure 7. Distribution of the "photometric members" of the clus- 
ter at z = 1. The area of the circles are proportional to the Ks- 
band brightness of the galaxies. The physical size of the radius of 
the large circle is 0.5 Mpc. The pixel scale is 0'.' 147 pix~ 1 . North 
is up, East is to the left. 



relations at the cluster redshift (z = 1 ) for passively evolv- 
ing g alaxies formed at Z/=2, 3 and 5 llKodama fc Arimotol 
1997). This model assumes that the red sequence of early 
type galaxies in the colour-magnitude diagram is a pure 
metalicity sequence as a function of galaxy magnitude, and 
calibrate its zeropoint using the red sequence of early type 
galaxies in the Coma cluster. 

Most of the photometric members are concentrated in 
the "inner region" around the lensing galaxy, while the non 
members are distribubuted more homogeneously across the 
field. The photometric members span a "red sequence" in 
the colour-magnitude diagrams with a considerable scatter. 
The slope and zeropoint of the red sequence are hard to 
define due to the large scatter but are roughly consistent 
with predictions of the simple models. 

The large observed scatter is likely to be caused by vari- 
ations in the dust content, the morphology and/or the age 
of the stellar populations of the galaxies. The red sequence 
is expected to increase its scatter and eventually fall apart 
as the redshift approaches the formation redshift of the stars 
in the galaxies. The properties of the red sequence are con- 
sistent with the picture of a young cluster in the process of 
formation, in which the cluster galaxies habour young stel- 
lar populations and have not yet evoloved in to the old red 
galaxies found in the cores of relaxed clusters. 



Figure 8. J-Ks versus Ks colour magnitude relation. Squares 
mark spectroscopic members. Large symbols are photometric 
members, while small symbols are non members. Galaxies within 
0.5 Mpc of the lensing galaxy are represented by filled symbols, 
while galaxies further away than 0.5 Mpc are represented by open 
symbols. The three lines are the predicted colour-magnitude rela- 
tions at the cluster redshi ft (z = 1) for passively evol ving galaxies 
formed at «/=2, 3 and 5 jKodama &: Arimotdll997l) . 




zf=5 
zf=3 
zf = 2 



20 
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Figure 9. R-Ks versus Ks colour-magnitude diagram for galax- 
ies detected in the R band. Symbols are as in Fig. [H] 



5.4 Cluster galaxy luminosity function 

We can now proceed to derive the Ks-hanA luminosity func- 
tion of the cluster galaxies. Since the T^s-band observations 
are complete to Ks = 22 and we know the photometric red- 
shift distribution of all galaxies in the field, we can do this 
without having to make uncertain statistical corrections to 
account for foreground and background contamination. To 



take full advantage of the data, we use a maximum likelihood 
technique applied directly to the luminosity distribution of 
the cluster galaxies, rather than doing a "least squares" fit to 
a binned representation. Details of the method are discussed 
in App^ In Fig. llUl we show the best fitting Schechter func- 
tion parameters a = -0.60lg;|| and K* = 18.90^;^. 

In Fig. llll we plot the derived Ks-band luminosity func- 
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Figure 10. Contour plot showing the constraints on the 
Schechter function parameters derived from the maximum like- 
lihood analysis. The cross marks the best fitting parameters, the 
curves represents the l-3a confidence levels 




Figure 11. ifs-band luminosity function of the cluster galaxies 
represented by the Schechter function with the a = — 0.60j^Qjg 
and K* = 18.90^Q jy. The filled symbols are a binned repre- 
sentation of the raw counts, while the open symbols have been 
corrected for incompleteness using the incompleteness function in 
Fig.0 



tion of the cluster galaxies. It is found to be well repre- 
sented by the Schechter function. The constraints on a are 
not strong, but it is important to leave it as a free parameter 
since the derived uncertainty on its value is coupled to the 
derived value of Ks* and its uncertainties. Also, since a has 
recen tly been shown to depend on wavelength jGoto et alJ 
2002) it is interesting to note that the value of a is in agree- 



Figure 12. Evolution of K* with redshift. The curves repre- 
sent the_£vohrtion_^red^ted_by stellar population synthesis mod- 
els iKodama &: Ari moto||l 997|). normalized to the Coma cluster 
which has Kt = 10.9 Ide^ropris et aflll99St) . Dashed lines are 
calculated in a Q m = 0.3, Ha = 0, ho = 0.6 cosmology while 
the full lines calculated in a Q m = 0.3, Oa = 0.7, ho = 0.7 cos- 
mology. Lines are labeled according to the formation redshift of 
their stellar populations. The lines labeled "No-evolution" only 
include K-corrections. The circles are the results of the present 
work. The filled circle is the result derived with a as a free pa- 
rameter, while the open circle is the result obtained with a fixed 
a = —0.9 . The diamonds and the square are resu lts from the 
literature ide Propris et alj[l999llNakata et al.ll200ll) 



ment with the value derived locally in the z*-band (which 
roughly corresponds to the same restframe wavelength as 
the Ks-band at z = 1) namely a z *,i oca i = —0.58. 

Earlier studies of the evolution of the luminosity func- 
tion with redshift have assumed the faint-end slope to be 
fixed at the value measured in the Coma cluster a = —0.9 
ide Propris et alll999HNakata et"al]|200lll . To allow a com- 
parison with these results we repeat our analysis with the 
faint-end slope fixed at a — —0.9 and derive a la con- 
straint on the characteristic magnitude of K* = 18.601q 24- 
In Fig. 1121 we compare the derived value of K* with values 
derived at other redshifts, and with the predictions of pas- 
sively evolving st ellar population models wit h different for- 
mation redshifts iKodama fc Arimotolll99^7l) . Dashed lines 
are calculated in a f2 m = 0.3, Oa = 0, ho = 0.6 cosmol- 
ogy while full lines are calculated in a fi m = 0.3, 57a = 0.7, 
ho = 0.7 cosmology. The filled circle is the value derived with 
a left as a free parameter. The open circle is the value de- 
rived with a fixed a = —0.9. We have included this result in 
the plot to make a more direct comparison with results from 
the literature. As expected from Fig. 1101 fixing a results in 
a brighter best fitting K* with smaller erro rbars, bringing it 
into p erfec t agreement with th e results of Ide Propris et ail 
( 1999) and Nak ata et alJ (2001) and consistency with what is 
expected for a passively evolving stellar populations formed 
at high redshift Zf > 2. 

To investigate the effects of field galaxy polution we re- 
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did the analysis on a subsample of the photometric member 
galaxies with photmetric redshifts in the range 0.9 < z p hot < 
1.1. This limits the effects of field galaxy polution but also 
excludes several known (bright) cluster galaxies from the 
sample. The best fitting values obtained in this way are 



-0-51±o.28 and K * 



19 - 10 ±o'.25- The favored it.,* is 



a bit fainter, due to the exclusion of some of the brightetest 
galaxies from the sample, but apart from that the main ef- 
fect is an increase of the size of the error bars. With the 
faint-end slope fixed at a — —0.9, the best fitting char- 
acteristic magnitude hardly changes. The derived value is 
it* — 18.551q' 27 . 

From the present analysis we conclude that the effects 
of field galaxy polution on the derived luminosity function 
paramenters are marginal, but that fixing the faint-end slope 
of the cluster galaxy luminosity function directly affects the 
derived characteristic magnitude. We find no evidence for 
a steepening or shift to fainter magnitudes of the bright 
end of the (cluster galaxy) luminosity function at 2 = 1, 
when compared to the local (z — 0) one. Nor do we find 
a significant change of its faint-end slope to support the 
hierachical formation scenario. 



5.5 Mass-to-light ratio 

The restframe it s -band mass-to-light ratio of the cluster 
around MG2016+112 can be estimated by comparing the 
total cluster mass, e stimated from the v elocity dispersion of 
6 cluster members (ISoucail et all 1200 ll) . with the total K- 
corrected its-band flux of the cluster galaxies. We calculate 
the absolute K-corrected its-band magnitude of the cluster 
galaxies to be 



M K s = m Ks — 51ogdi - 25 — k Ks (z = 1), 



(1) 



where d,L is the luminosity distance in Mpc calculated in 
a Q m — 1, Qa = 0, ho = 0.5 (for easier compariason to 
previous results), and kics{z = 1) = —0.352 is the k correc- 
tion in th e its band for an object at redshift z = 1 from 
IPoggiantl jl997f) (we apply this term to be able to com- 
pare with the absolute its-band magnitude of the sun. We 
assume M v q = 4.8 and (V — Ks)q = 1.45). 

The derived its-band mass-to-light ratio of the clus- 
ter is found to be M/L Ks = 27± f $ h 50 (M/L Ks )q. This 
is in agreement with the upper limi t on the total mass 
of the cluster derived from X-rays JChartas et alJ l200ll) 
which when combined with our estimate of the total its- 
band luminosity of the cluster galaxies translates into 
the following upper limit on its-band mass-to-light ratio 
M/Lkb < 43/i5o (M/ Lks)q, and in agreement with the 
its-band mass-to-light ratio derived lo cally in the Com a 
cluster M/Lks = 38±23/i 5 o (M/L Ks ) Q iRines et all200ll) . 



6 SUMMARY 

We have obtained deep NIR and optical imaging of the field 
around the wide separation gravitational lens MG2016+112. 
A photometric redshift analysis of objects in the field, reveal 
69 galaxies with photometric redshifts consistent with being 
in a cluster at the redshift of the lensing galaxy 2=1. 

The its-band luminosity function of the cluster galax- 
ies is well represented by the Schechter function. The best 



fitting faint-end slope a = — O.BOig'^ is consistent with 
what is measured at the same restframe wavelength (the 
z* band) locally, providing no evidence of evolution. In the 
bright end no steepening or shift to fainter magnitudes of 
the luminosity function is found. The best fitting charac- 
teristic magnitude it* = 18.90lo'g 7 is consistent with what 
is expected for a passively evolving population of galaxies 
formed at high redshift Zf > 2. From the its-band flux 
of the cluster ga l axies and the total cluster mass estimates 
of ISoucail et alJ i200ll) a its-band mass-to-light ratio of 
M/Lks = 27^17/150 [M/ Lks)q is derived, consistent with 
the upper limit deri ved from Chandra X-ray observations 
JChartas et alJl200lh and the va lue measured locally in the 
Coma cluster 1 Rines et al.ll200lh . 

The projected spatial distribution of the cluster galaxies 
is filamentary-like rather than centrally concentrated around 
the lensing galaxy, and show no apparent luminosity segre- 
gation. A handful of the cluster galaxies shows evidence of 
merging and interaction. The cluster galaxies span a red se- 
quence with a considerable scatter in the colour-magnitude 
diagrams, suggesting that they contain young stellar popu- 
lations in addition to the old populations of main sequence 
stars that dominate the its-band luminosity function. This 
is in agreement with spectroscopic observations which show 
that 5 out of 6 galaxies at the redshift of the lensing galaxy 
has emission lines. 

The evidence presented here suggest that a young clus- 
ter of galaxies is in the process of assembling around the 
massive lensing galaxy MG2016+112. The cluster is mass- 
selected (from the precence of the wide-seperation gravita- 
tionally lensed QSO) and is therefore likely to be more typi- 
cal of the z ~ 1 cluster population than clusters selected us- 
ing X-ray or optical techniques that rely on the presence of a 
relaxed intra cluster medium and / or a central concentration 
of old bright red galaxies for the clusters to be detected. 



APPENDIX A: DERIVING THE CLUSTER 
GALAXY LUMINOSITY FUNCTION USING 
THE MAXIMUM LIKELIHOOD TECHNIQUE 

To take full advantage of the data, we apply th e maximum 
likelihood technique of ISchechter fc Pressl il976f) directly to 
the luminosity distribution of the cluster galaxies, rather 
than doing a "least squares" fit to a binned representation. 
In this way we include the information that in many mag- 
nitudes intervals no galaxies are found, and we do not make 
the erroneous assumption of the \ 2 method that the under- 
lying distribution is gaussian (it is in fact poissonian, either 
a galaxy is there or it is not). This in turn leads to results 
with more realistic errorbars. 

In the magnitude interval between m and m + dm, the 
expected number of galaxies p(m)dm is taken to be 



p(m)dm = f(m)n* I0°' 4 ' m m ' 



a + l 



dm, (Al) 



where m* is the characteristic magnitude of the cluster 
galaxies, a is the slope of the faint-end of the cluster galaxy 
luminosity function, n* is a normalization constant and 
f(m) is the completeness function of the "photometric mem- 
ber" sample (shown in Fig.QjJ evaluated in m. 

Poisson statistics gives the probability that a cluster of 
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galaxies characterized by n* , m* and a will yield the data 
set actually observed (N; 772l, 7722 5 — ,772jV , l^tlim) 



P = exp 



p(m)dm 



[p(m;)dm] 



(A2) 



where mi» is the magnitude limit (to which the dataset is 
complete, in this case mu m = Ku m — 22). The parameters 
a and m* which are "most likely" to represent the observed 
data set are tho se that min imizes the C statistics (the like- 
lihood function dCashlll979l) 1 



C = -21n(P). 



(A3) 



The minimum value of the C statistics Cmm is found by 
varying a and m* . For each combination of on and m*,i 
we calculate the term AC = C(ai,m„ t i) — C m j n , which i s 
distributed as \ 2 with 2 degrees of freedom l)Cashlll979h . 
From AC the l-3u confidence levels can be derived in a 
straight forward manner as these (as for the \ 2 distribution) 
are found where AC = 2.29575, 6.18008, 11.8291. 
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